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Fiber  composite  materials  have  very  high  specific 
strength  and  stiffness.  This  structural  efficiency  enables 
composite  materials  CO  be  used  in  a wide  variety  of 
applications.  However,  laminated  composites  have  also  some 
disadvantages  such  as  poor  inter- laminar  strength,  low 
iagsact  resistance,  and  poor  deXamlnation  resistance.  During 
the  manufacture  of  fiber-reinforced  composite  laminates, 
imperfections  can  cause  initial  delaminatlon.  and  the 
impact  on  laminated  composites  by  foreign  objects  during 
service  can  cause  delaminatlon  propagation.  Delamination  is 
one  of  the  most  dominant  forms  of  damage  due  to  lack  of 
reinforcement  in  the  thickness  direction.  An  effective 


method  to  provide  this  reinforcement  is  to  use  through-the- 
thickness  stitching. 

Hhen  the  delaminated  composite  is  under  cor^ression, 
it  may  undergo  post-buckling,  nhen  the  energy  release  rate 
is  greater  than  critical  energy  release  rate,  the 
delamination  will  propagate  in  an  unstable  manner  and 
finally  cause  catastrophic  failure  of  the  composite.  In 
this  dissertation  analytical  and  finite  element  methods  are 
proposed  to  investigate  the  effect  of  through-the-thlckness 
stitching  on  the  post-buckling  of  a delaminated  composite. 
In  particular  the  studies  focus  on  the  effect  of  stitching 
on  buckling  load,  load-displacement  behavior,  and  energy 
release  rate. 

When  a delaminated  composite  is  subject  to  a low- 
velocity  intact,  the  crack  may  propagate  and  finally  cause 
catastrophic  failure.  In  this  dissertation  an  analytical 
model  is  proposed  to  understand  the  effect  of  stitching  on 
improving  impact  damage  resistance.  In  particular  the 
studies  focus  on  the  effects  of  stitching  when  the  impact 
energy,  delamination  length  and  location,  and  the  fracture 
toughness  are  varied. 


CHAPTER  1 

INTRODDCTION  AMD  IITERATUBE  REVIEM 
1.1  Bacliqtound 

Fiber  composite  materials  have  very  high  strength  and 
stiffness  hut  low  density.  This  structural  efficiency 
enables  these  composite  materials  to  be  used  in  a wide 
variety  of  applications,  especially  in  aerospace 
structures.  However,  advanced  composites  also  have  some 
disadvantages  such  as  poor  inter-laminar  strength,  low 
damage  resistance,  low  damage  tolerance,  and  poor 
delamination  resistance.  Since  most  of  the  commonly  used 
advanced  composites  are  made  by  laying  up  fibrous 
reinforcement  in  a surrounding  matrix  and  the  mechanical 
properties  between  the  plies  are  dominated  by  a relatively 
weah  matrix,  this  weakness  causes  low  damage  resistance  and 
low  damage  tolerance  and  finally  results  in  low 
eompresslon-af cer-lnpact  strength. 

A delamination  is  an  interface  crack  or  a debond 
between  two  adjacent  plies.  It  is  one  of  the  critical 
failure  modes  in  composite  materials.  During  the 
manufacture  of  fiber-reinforced  composite  laminates, 


imperfections  such  es  air  entrapment  or  resin  starvation 
may  cause  initial  delanination  in  the  con^site.  Impact  on 
laminated  composites  by  a foreign  body  during  service  may 
also  cause  delamination  while  the  surface  of  the  composite 
structure  remains  undamaged  to  visual  inspection.  This  kind 
of  delasdnation  can  randomly  occur  somewhere  in  the 
structure.  No  matter  where  it  occurs  it  may  decrease  the 
overall  stiffness  and  load-carrying  capacity  of  the 
composite  structure.  The  growth  of  delamination  cracks 
under  the  subsequent  application  of  external  loads  leads  to 
the  rapid  deterioration  of  mechanical  properties  and  may 
cause  catastrophic  failure  of  the  composite  structure.  One 
of  the  most  dominant  forms  of  damage  is  delamination  due  to 
lack  of  reinforcement  In  the  thickness  direction.  There  are 
some  approaches  that  can  be  adopted  in  order  to  improve  the 
strength  in  the  thickness  direction.  Using  fully  Integrated 
3D  composites  such  as  weaves  and  braids  is  one  of  the 
methods;  however,  due  to  their  complexity,  limited  shape- 
ability.  and  processability,  their  applications  have  been 
limited.  Inserting  a pin  in  the  thickness  direction  can 
also  improve  reinforcement  in  the  thickness  direction,  but 
currently  it  is  difficult  to  use  this  technique  in  large 
structures.  Another  method  for  providing  a reinforcement  in 
the  thickness  direction  is  to  use  stitching  through  the 


thicluiess  of  the  lasinate.  Only  small  amounts  of  out-of- 
plane  reinforcement  are  needed  in  order  to  significantly 
change  the  mechanical  properties  of  the  structure.  A 
stitched  laminate  can  be  formed  by  using  a sewing  machine 
to  stitch  fiber  yarns  into  a composite  preform  before 
curing.  Having  a stitch  in  a laminate  does  not  greatly 
alter  the  original  laminated  structure.  Further,  stitching 
can  be  done  very  fast  in  automated  machines,  and  is 
suitable  for  large  structures  such  as  aircraft  wings. 

Under  compression  the  delaminated  part  of  composite 
laminate  may  undergo  post-bucbiing  behavior.  A delamination 
may  propagate  when  a foreign  body  drops  on  a delaminated 
composite.  Many  researchers  have  studied  these  subjects  by 
using  analytical  and  numerical  methods.  Some  studies  have 
investigated  the  effect  of  stitching  on  mechanical  and 
Impact  properties,  such  as  inter-laminar  strength  and 
fracture  toughness. 

A literature  review  ie  provided  in  the  next  section. 

1.2  Literature  Review 

The  literature  review  is  divided  into  three  partsr  (a) 
post-buchling  behavior  of  delaminated  composites;  (b) 
impact  response  of  composite  structures;  and  (c)  through- 
the-thicbness  stitching  of  composite  laminates. 


1.2.1  Research  on  Post-bucKling  of  Delaminated  Composites 

In  nany  caseS/  when  delaminated  coinposites  are  under 
external  compressive  loads  the  delaminated  structure  will 
have  Che  capability  of  carrying  loads  beyond  their  buckling 
loads  and  may  fall  in  the  post-buckling  regime.  Therefore 
bifurcation  analysis  may  not  be  adequate  in  describing  the 
whole  damage  process  of  delaminated  composites  up  to  final 
failure;  a post-buckling  analysis  may  be  required. 

The  Initial  development  and  growth  of  delanlnation  in 
compressed  composites  are  affected  by  various  geometrical 
parameters,  loading  conditions,  material  properties  and  so 
on.  The  different  combinations  of  these  parameters  can 
result  in  different  types  of  buckling  and  delaminatlon 
growth  behavior.  Host  studies  were  based  on  the  static 
post-buckling  solution  of  homogeneous  or  laminated  plate 
theory.  The  energy  release  rate  (C)  is  considered  as  the 
critical  parameter  and  different  methods  were  used  to 
calculate  the  energy  release  rate.  An  assumption  that  uses 
Griffith-type  fracture  criterion  with  a specific  constant 
critical  fracture  energy  release  rate  (G^l  was  adopted  in 
these  studies.  When  the  energy  release  rate  S is  greater 
than  critical  fracture  energy  release  rate  G«  the 
delamination  in  the 


composite  may  propagate. 


Slciple  analytlcsl  and  finite  element  models  seem  to  be 
popular  in  analysing  one-dimensional  delaminations . Chai  et 
al.  [11  used  an  analytical  method  to  obtain  energy  release 
rate  by  differentiating  the  strain  energy  with  respect  to 
delamination  length,  fin  and  Mang  |2]  used  the  J-lntegral 
method  to  obtain  an  expression  for  energy  release  rate  in 
the  post-buckled  delaminated  plate.  Hhitcomb  [31  adopted 
the  crack  closure  technique  and  analyzed  the  effects  of 
various  parameters  on  energy  release  rate  in  a beam-like 
structure  using  geonetrically  nonlinear  finite  elements. 

Sheinman  and  Softer  [4]  investigated  the  effects  of 
bending-extension  coupling  and  initial  imperfection  on  the 
post-buckling  behavior  of  a composite  beam  with  various 
delamination  geometries.  It  was  found  that  the  coupling 
effect  significantly  reduced  the  buckling  load  and 
increased  the  post-buckling  deformations,  and  the  global 
post-buckling  deformation  was  shown  to  be  very  sensitive  to 
initial  in^erfections. 

Kardomateas  and  Schmueser  [S]  proposed  a one- 
dimensional  beam-plate  model  that  accounts  for  the 
transverse  shear  effect.  Analytical  solutions  for  the 
critical  instability  load  and  the  post-buckling  deflsctlons 
were  obtained  with  aid  of  the  perturbation  technique. 
Kardomateas  [6]  also  modeled  the  post-buckliog  behavior  of 


a laminate  with  a thin  delaminated  aublaminate  through  a 
procedure  that  is  based  on  large  deflection  theory  of  the 
delaminated  layer. 

Chen  ['7-61  proposed  a shear  deformation  theory  and 
obtained  energy  release  rate  by  a variational  energy 
principle.  He  found  that  the  energy  release  rate  is  larger 
with  transverse  shear  effect  Included  for  the  same  applied 
load.  Therefore  the  initiation  of  delamination  growth  will 
occur  at  a lower  applied  load  than  that  evaluated  with  the 
classical  lamination  theory.  He  also  found  that  the 
magnitude  of  the  transverse  shear  effects  depends  on  the 
delamination  location  and  sire. 

Davidson  and  Krafchak  |9],  Davidson  (ICI,  Davidson  and 
Krafchah  [11]  used  crack  tip  element  analysis  to  determine 
total  energy  release  rate  end  Individual  mode  1 and  II 
energy  release  rate  in  different  loading  conditions  or 
different  delamination  cases. 

Kyoung  and  Kim  [12|  presented  an  energy  method  to 
obtain  an  analytical  solution  for  determining  the  buckling 
load  and  the  growth  of  delaminated  beam-plate  structures. 

In  their  study  the  delamination  was  arbitrary  located.  They 
elso  found  that  the  shear  effects  decrease  the  buckling 
loads  and  increase  energy  release  rate.  The  results  from 
their  analysis  show  that  the  buckling  load  of  a center 


located  delajulnation  is  higher  than  an  off-center 
de lamination. 

Lee  et  al.  [13|  proposed  a finite  element  method  based 
on  a layer-wise  laminated  composite  plate  theory  for 
solving  the  post-buckling  problem.  In  their  study  a contact 
algorithm,  which  overcame  the  physically  unallowable 
overlapping  between  delaminated  surfaces,  was  used. 

There  are  also  several  works  concerning  two  dimension 
post-buckling  analysis.  Yin  and  Fei  [14]  presented  an 
elastic  post-buckling  analysis  of  a delaminated  circular 
plate  under  axis-symmetric  compression  along  Its  clamped 
boundary.  They  found  that  certain  features  of  the  post- 
buckling  behavior  are  quelitatlvely  similar  to  the  buckling 
of  an  axially  loaded  beam  plate  containing  a one- 
dimensional delamination.  Whitcomb  and  Shivakumar  [151  used 
virtual  crack  closure  technique  to  calculate  Che  total 
energy  release  rate  In  a locally  post-buckled  laminate  with 
embedded  delamination.  They  found  that  there  is  a large 
variation  of  energy  release  rate  along  the  delamination 
front  for  square  and  rectangular  plates,  and  delamination 
growth  depends  on  current  delamination  aspect  ratio,  the 
strain  level,  and  the  size  of  the  delamination. 

Suemasu  [16]  Investigated  an  effect  of  multiple 
delaminations  of  composite  panels  on  the  buckling  loads  by 


using  an  analytical  approach  and  verified  the  loodel  using 
experiments  and  finite  element  methods.  He  pointed  out  that 
the  buckling  load  reduces  significantly  due  Co  the 
existence  of  multiple  delarainations . 

Naganarayana  and  Atluri  [11-191  used  a multi-plate 
model,  in  conjunction  with  a 3-noded  quasi-comforming  shell 
element,  Co  model  the  delaminated  plates.  The  J-integral 
technique  was  used  for  delamination  growth  prediction  in 
term  of  point-wise  energy  release  rate  distribution  along 
the  delamination  edge.  The  effects  of  structural  parameters 
such  as  delamination  thickness,  size  and  shape  on  the  post- 
buckling  behavior  and  on  Che  delamination  growth  were 
examined. 

Sben  and  Hiliiams  (19]  studied  a post-buckling 
analysis  of  laminated  plates  subjected  to  combined  axial 
compress  and  uniform  temperature  loading.  The  analysis  used 
a perturbation  technique  to  determine  buckling  loads  and 
post-buckling  behaviors  in  anti-symmetrically  angle-ply  and 
symmetrically  cross-ply  laminated  plates.  Klug  et  al  [20] 
adopted  Mindlin  plate  finite  element  Co  perform  the  post- 
buckling  analysis  and  to  compute  energy  release  rate  at  the 
delamlnatlon  front  with  aid  of  the  crack  closure  method  in 
composite  laminates  which  have  elliptic  delamination  shape. 
The  results  show  Chat  this  method  is  efficient  and  accurate 


compared  to  three-dimensional  finite  method  for  calculation 
of  the  energy  release  rate. 

1.2.2  Research  on  Low-Velocity  Impact  Reaponse  and  Damage 
of  Delaminated  Composite  Materials 

Low-velocity  impact  response  and  damage  caused  by  low- 
velocity  projectile  impact  has  received  attention  from  some 
researchers.  When  the  low-velocity  foreign  body  drops  on  a 
structure  such  as  the  surface  of  an  airplane  wing,  it  may 
cause  damage  such  as  a delamination  which  is  invisible  on 
the  surface  and  difficult  to  detect.  Therefore 
understanding  the  low-velocity  impact  response  of  a 
delaminated  composite  is  very  important,  and  a method  is 
needed  to  predict  the  damage  due  to  low-velocity  impact. 
Graves  and  Kantr  1211  used  the  maximum  shear  stress  failure 
criterion  to  estimate  the  damage  area.  Clark  |22|  provided 
a model  only  for  qualitative  prediction  of  the  delamination 
sire.  Qrady  and  Sun  1231.  and  Grady  and  Depaola  (24] 
provided  an  estimation  of  the  delamination  growth  due  to 
impact.  Choi  et  al.l25|  adopted  a dynamic  finite  element 
method  associated  with  failure  analysis  to  predict 
threshold  of  impact  of  damage  and  initiation  of 
delamination.  Sankar  and  Hu  126]  and  Hu  [27]  used  a finite 
element  method  to  analyse  dynamic  delamination  growth  in  a 


composite  beam.  The  delaminated  beams  were  modeled  as  two 
offset  beams  and  spring  elements  were  used  to  connect  the 
beams  in  the  uncrac)ced  portion.  The  dynamic  energy  release 
rate  was  computed  from  the  strain  energy  in  Che  crack'tlp 
element.  Razi  and  Kobayashi  |2B]  performed  finite  element 
analysis  and  experiments  Co  study  damage  growth  and 
distribution  In  cross-ply  laminate  beams  and  plates.  A 
guasi-static  finite  element  analysis  coupled  with  a Gjic 
criterion  was  proposed  for  the  delamination  growth/  and  Qui 
Gut  criterion  was  used  to  arrest  crack  propagation.  Sun  and 
Yih  [291  investigated  the  quasi-static  characteristics  of 
impact  responses  and  the  impact-induced  delamination  growth 
of  composite  laminates  subjected  Co  low-velocity  impact. 

They  used  a spring-mass  model  to  predict  contact  force 
history  and  the  peak  force  was  associated  with  the  energy 
release  rate.  Thermal  loads  were  included  in  computing 
energy  release  rate.  Wang  and  Vu-Khanh  (30]  also  used  the 
finite  element  method  to  model  carbon  fiber/PEEK  cross-ply 
damaged  laminates.  They  found  that  a delamination  occurred 
in  Che  form  of  a Mode  ll-dominated  unstable  crack  growth 
and  subsequent  arrest. 

Kwon  and  Sankar  131-32]  performed  experiments  to  study 
the  low-velocity  impact  damage  of  quasi-isotropic  and 
cross-ply  graphite/ epoxy  composite  laminates.  They  found 


that  tha  impact  force  history  and  delamination  radius  could 
be  predicted  using  the  results  from  a static  indentation- 
flexure  test.  Sanhar  t33]  reviewed  the  necessary  concepts 
for  understanding  low-velocity  impact  response  and  damage 
in  fiber  composite  materials  and  gave  several  algorithms 
for  predicting  the  impact  force  history  to  different 
degrees  of  approximations.  Sankar  [331  presented  a semi- 
empirical  method  for  predicting  impact  damage  in 
composites. 

1.2.3  Research  Related  to  Delaminated  Composite  Haterials 


In  transversely  stitched  composites  the  delamination 
is  not  completely  separated  but  is  held  together  by  the 
stitches.  Therefore  the  propagation  of  cracks  can  be 
prevented  in  stitched  composite  materials.  Several  studies 
have  investigated  the  effect  of  stitching  on  mechanical  and 
impact  propertiss,  failure  loads,  inter-laminar  strength, 
and  fracture  toughness. 

In  general  there  are  three  types  of  stitches  as  shown 
in  Figure  1.1;  modified  lock  stitch,  chain  stitch,  and 
standard  lock  stitch.  Some  parameters  of  stitches  are: 
Stitch  density,  stitch  yarn,  and  yarn  material  density.  We 
define  the  stitch  density  in  a composite  laminate  by  the 
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number  of  stitches  per  square  inch  and  represent  this 
density  by  the  stitching  pattern  as:  (number  of  stitches 
per  inch)  X (spacing  between  two  stitch  lines),  for 
example,  4X1/4  means  a stitch  density  of  16  where  there  are 
4 stitches  per  inch  and  distance  between  two  adjacent 
stitch  rows  is  1/4  in.  Stitch  yarn  is  the  bobbin  yarn  used 
in  the  stitching  process.  The  material  density  of  the 
stitch  yarn  is  expressed  as  denier.  One  denier  is  the  mass 
in  grams  of  9000  meters  of  yarn. 


rigure  1.1  Types  of  locJc  stitches 
(Courtesy,  Sharma  and  Sanlcar  (461) 

Sawyer  (34]  used  an  experimental  method  to  determine 
the  effect  of  stitching  on  the  static  load  of  bonded 
composite  single  lap  joint.  Up  to  334  improvement  in  static 
failure  load  compared  with  unstitched  results  was  obtained 
by  a single  row  of  stitches  near  the  end  of  the  overlap. 


Additional  row  of  stitching  or  different  stitch  spacing  has 
little  effect  on  static  joint  failure  load. 

Su  [351  Improved  the  delamination  resistance  in 
composite  laminetes  by  the  use  of  thermoplastic  matrix 
resin  and  stitching.  There  was  a 20-30i  further  improvement 
of  delaminetion  resistance  by  using  stitches.  Dexter  and 
Funk  (361  conducted  an  experimental  investigation  on  impact 
resistance  and  inter-laminar  fracture  toughness  of  stitched 
guasi-isotropic  graphite/epoxy  composites.  The  stitch 
materials  were  polyester  or  Kevlar""  yarns  with  several 
parameters.  They  found  that  a significant  drop  in  damage 
areas  of  stitched  laminates  compared  to  unstitched 
laminates  for  the  same  impact  energy.  The  results  shown 
that  inter-laminar  fracture  toughness  Gis  of  stitched 
laminated  with  the  most  effective  Kevlar  yarn  was  30  times 
higher  than  unstitched  laminates. 

Pelstring  and  Madan  [37]  proposed  a semi-empirical 
formula  to  relate  a damage  tolerance  of  a composite 
laminate  to  stitching  parameters.  The  results  showed  that 
Mode  I critical  energy  release  rate  of  a stitched  laminate 
was  15  times  higher  than  an  unstitched  laminated. 

Dransfield  et  al.  (38)  reviewed  the  effect  of  stitching  on 
improving  some  mechanical  properties  such  as  compress- 
®fter-impact  strength  and  delamination  resistance  of 


composite  materials.  The  advantages  and  disadvantages  of 


adopting  a stitch  In  a composite  structure  and  the 
necessary  techniques  of  manufacture  were  examined.  Shu  and 
Mai  [39-40]  used  analytical  approaches  to  investigate  Che 
effect  of  stitching  on  buckling  load  and  energy  release 
rate  of  thin  film  structure.  They  assumed  that  Che  stitches 
followed  a Winkler  elastic  foundation  type  of  stress- 
separation  relation.  The  results  shown  that  the  strength  of 
composite  laminates  under  edgewise  compression  can  be 
significantly  increased.  Jain  and  Mai  [411  proposed  two 
analytical  models  to  study  the  effect  of  stitch  parameters 
such  as  stitch  density,  stitch  thread  diameter,  and 
matrix/stitch  inCerfacial  shear  stress  etc.  on  improving 
Mode  I delamination  toughness.  They  found  that  a large 
matrix/stitch  intertacial  shear  stress,  high  stitch  density 
and  a small  stitch  thread  diameter  were  desirable  in  order 
to  maxlmlke  the  delamination  growth  resistance.  Jain  and 
Hai  [42]  also  used  an  analytical  method  Co  study  the  effect 
of  stitching  parameters  on  Mode  II  energy  release  rate  in  a 
laminated  composite  and  found  that  through-thickness 
stitching  can  improve  Mode  ll  fracture  toughness  also.  Chen 
and  Sun  [43}  used  a finite  element  method  to  propose  an 
effective  energy  release  rate  for  measuring  Mode  I fracture 
toughness  of  a stitched  composite  laminate.  Kang  and  Lee 


Iflfll  did  an  experiment  to  study  the  effect  of  stitch  on 
mechanical  and  impact  properties  of  woven  laminated 
composites.  The  mechanical  properties  were  improved  at  an 
optimum  stitch  density  compared  with  unstitched  woven 
laminated  compoeites.  The  damage  area  caused  by  repeated 
impact  was  tar  smaller  in  stitched  woven  composites  than 
that  in  unstitched  composites.  Wu  and  Liau  (45)  also  uaed 
an  experimental  method  to  investigate  the  behavior  of 
stitched  laminates  under  low  velocity  impact.  They  also 
found  that  stitching  could  significantly  reduce  the 
delamination  area.  With  increasing  the  stitch  density  the 
damage  mechanism  changed  from  delamination  in  unstitched 
laminates  to  the  plastic-hinge  type  of  local  deformation  at 
the  impact  location  in  stitched  laminates.  Sharioa  and 
Sankar  146)  performed  experiments  to  study  the  effects  of 
stitching  on  CAI  (Compression-After-lmpact)  strength, 
fracture  toughness  in  Mode  I and  Mode  II.  They  found  that 
CAI  strength,  Fracture  toughness  of  Mode  I and  Mode  II  can 
be  improved  significantly  by  stitching.  Sankar  and  Sonlk 
(47)  used  Finite  Element  and  analytical  models  to  study  the 
effect  of  stitching  on  Mode  II  fracture  toughness. 

Dickinson  C4B]  adopted  a micro-mechanics  approach  to 
characterise  the  effects  of  stitching  on  the  elastic 
constants,  in-plane  strength  and  inter-laminar  mechanical 


response  of  undemaged  composite  laminates.  Sankar  and 
Dhamapuri  149]  presented  an  analytical  mode  to  describe 
delamlnation  growth  in  laminated  DCB  {Double  Cantilever 
Beam)  specimen  with  stitches  and  proposed  a simple  method 
to  compute  Che  energy  release  rate  doe  to  delamlnation  in 
the  stitched  specimen.  Jain  and  Mai  [50]  studied  Che  effect 
of  reinforced  Cabs  on  the  Mode  I fracture  toughness  of 
stitched  composites.  They  indicated  that  reinforcing 
aluminum  tabs  along  the  length  of  stitched  DCB  specimens 
can  alter  the  failure  mechanism  of  the  stitch  threads  and 
hence  erroneous  values  for  Mode  I fracture  toughness. 
Glaessgen  eC  al  [51]  used  finite  element  method  In 
conjunction  with  crack  closure  technigue  to  study  the 
effect  of  stitching  on  Mode  I and  Mode  II  energy  release 
rate  for  various  debond  configurations. 

1-3  Objectives  and  Scope  of  the  Present  Study 

From  the  above  sections  we  can  realize  that  one  of  the 
most  dominant  damage  in  composites  is  delamlnation  growth 
due  to  lack  of  reinforcement  in  the  thickness  direction. 
Stitching  is  found  to  be  one  of  the  effective  reinforcing 
methods  in  the  thickness  direction  of  laminated  composites. 
However  it  Is  not  clearly  understood  how  stitching  can 
delay  the  delamlnation  propagation  of  composite  laminates 
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subjected  to  compressive  loads  or  to  low-velocity  impact. 
The  purpose  of  the  present  research  is  to  use  analytical 
end  finite  element  methods  to  study  the  effect  of  stitching 
on  posC-buclcling  and  low-velocity  impact  response  of 
composite  structures.  In  particular  the  study  will  focus  on 
the  effect  of  stitching  on  buckling  load,  energy  release 
rate,  and  Impact  damage  resistance. 

The  organization  of  this  dissertation  is  as  follows: 

In  Chapter  2 an  analytical  solution  tor  post-buckling  of  a 
delaminated  composite  beam  with  stitching  will  be  derived, 
In  Chapter  3 a finite  element  analysis  will  be  given 
related  to  post-buckling  of  a delaminated  composite  beam 
with  stitching.  In  Chapter  4 an  analytical  model  tor  impact 
response  of  stitched  composite  beam  is  presented. 
Conclusions  drawn  from  this  study  and  discussions  for 
future  research  are  presented  In  chapter  5. 


CHAPTER  2 

AN  ANALYTICAL  SOLUTION  FOR  POST-BUCKLING  OF  A DELAMINATED 
CCWPOSITE  BEAM  WITH  STITCHING 


In  this  chapter  an  analytical  solution  for  post- 
buckling  behavior  of  a delaminated  composite  beam  with 
stitching  is  derived  baaed  on  Rayleigh-Bitz  method.  The 
effect  of  stitching  on  buckling  load  and  energy  release 
rate  is  discussed.  Non-dimensional  parameters  for  buckling 
load,  energy  release  rate  (G),  crack  length  (a),  and  stitch 
stiffness  (K)  are  derived.  Although  finite  element  methods 


can  model  a stitched  composite  more  realistically, 
analytical  methods  have  several  advantages.  Analytical 
methods  provide  a closed-form  solution  to  all  the  desired 


quantities  such 
buckling  loads. 


Ass^upptiotis 


When  we  model  Che  p 
deldmlnaCed  composite  be 
essumpCions  ace  made: 


iC-buckllng  behavior  of  a 
1 with  stitching  the  following 


..A  stitch  is  inodeled  as  a linear  cohesive  spring.  The 
stiffness  of  the  spring  is  a function  of  stitch 
density,  young's  modulus  £ of  the  stitch,  diameter 
of  the  stitch  yarn,  and  thickness  of  Che  stitched 
composite  beam.  Stitches  are  smeared  averaged  on  the 
whole  length  of  Che  beam,  and  thus  the  smeared 
stitch  stiffness  should  account  for  Che  stitch 
spacing  also. 

!.  A delamination  is  assumed  to  be  pre-existing  before 
Che  compressive  loading  Is  applied  and  is  located  at 
the  center  of  Che  beam  both  in  length  direction  and 
in  thickness  direction.  Therefore  the  beam-llke 
structure  is  symmetric  (see  Figure  2.11. 

1.  The  deformation  shape  of  the  stitched  composite  beam 
under  post-buckling  behavior  ia  assumed  to  be  the 
same  as  in  an  unstitched  beam  under  post-buckling. 

. The  laminated  beam  is  approximated  as  a homogenous 
beam  with  an  equivalent  Young's  modulus. 


Analytical  Model 


In  this  section  we  first  derive  the  buckling  load  Per 
and  energy  release  rate  (G)  in  an  unstitched  composite  beam 
using  Che  Rayleigh-Ritz  method  and  Chen  obtain  the  buckling 
load  Per  and  energy  release  rate  (G)  for  a stitched 
composite  beam.  A single  crack  is  located  at  the  center  of 
the  structure. 

2.2.1  Delaminated  Composite  Beam  Without  Stitching 

A composite  beam-llke  structure  without  considering 
through-Che-Chickness  reinforcement  is  shown  In  Fig.  2.1. 


Figure  2.1  A delaminated  beam 

The  length  of  the  beam  is  2£,  S is  Che  width,  2h  is  the 
thickness  of  Che  beam,  2a  is  the  delaminaCion  length  and  P 
is  the  compressive  load. 


(2Z) 


Here  Uo  Is  the  raid-plane  displacecoent  in  the  x direction,  w 
is  the  displacement  of  top  sublamated  bean  In  the  a 
direction,  e,  is  the  normal  strain  of  top  sublaminated  beam 
in  X direction.  The  cross-sectional  area  of  the 
sublaminates  is  represented  by  A and  E is  the  Young's 
modulus.  Assume  the  forms  of  Ug  and  w as;  uo«cx, 
w=dsin'nxy2a,  where  c and  d are  unknown  coefficients  to  be 
determined. 

The  strain  energy  of  the  delaminated  portion  of  the 
beam  can  be  derived  as: 


(2  3) 


sublaolnate. 


The  potential  energy  V of  the  external  loads  is  given  by 


f'  = -(-/>)r(2a)  (2.4) 

Total  energy  FI: 

n = u+v 


Using  the  Raylelgh-Rit*  method: 


Substituting  tor  the  energy  terms  U and  V we  can  obtain  the 
fallowing  two  equations: 


(2  8) 


fs’r  ^ 3a-V’|  AFJn 
lie  8{2d)’J‘^  (2o)’ 


(2  9) 
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It  may  be  noted  that  the  above  equations  are  nonlinear 
and  lead  to  multiple  solutions.  Solving  these  two  equations 
we  obtain  the  coefficients  c and  d as  follows: 

One  passible  solution  set  is: 


(2.10) 


This  solution  describes  the  behavior  of  the  cracked 
portion  of  the  beam  where  the  compressive  load  is  less  than 
the  buckling  load  Per-  Since  d=0,  there  is  no  deformation  in 
1 direction.  The  relation  between  compressive  load  P and 
axial  displacement  ut  takes  a simple  form: 

Px 

“•*  2Fji  (2.12) 


The  end  shortening  A at  x=2a  is  given  by: 


The  other  possible  solution  set  is; 


(213) 


3P 

2AE 


Sx'l 

A(2o)' 


(214) 


25 

*P{.to)'  32/ 

KAn^  A 


(2  15) 


/’>P„ 


la^El 


(2  16) 


here  Per  is  the  critical  buckling  load. 

This  solution  set  describes  the  behavior  of  the 
cracked  portion  of  the  beam  where  the  compressive  load  is 
greater  than  the  buckling  load  P„.  There  is  a deformation 
in  the  2 direction.  The  relation  among  compressive  load  P 
and  axial  displacement  Uc  and  deformation  shape  w ace: 


(2,17) 


(2.18) 


here  d is  obtained  from  Equation  2.15  for  a given  P.  End- 
shortening  6,  at  x=2a  can  be  expressed  as  shown  as: 


^](2o 

2a?  f 


(2.19) 
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Based  on  Equation  2.13  and  2.19  the  relation  between  the 
compressive  load  P and  the  end  shortening  iwill  be  as 
shown  in  Figure  2.3. 


Figure  2.3  The  relation  between  compressive 
load  F and  end  shortening  6, 

In  Figure  2.3  S,r  is  end  shortening  at  critical  budding  load 
Per  and  is  end  shortening  at  an  arbitrary  compressive 

Based  on  Figure  2.3,  the  strain  energy  U,  of  the 
cracked  portion  can  be  written  as: 

3/>|'o  As‘EI’ 

‘^£4  Ar' 


(2  20) 
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He  can  extend  the  cesulta  to  the  entire  oonposite 
bean.  After  the  compressive  loads  P reaches  the  buckling 
load  Pei,  the  bean-like  structure  will  deform  as  shown  in 
Figure  2.4 , 


The  relation  between  the  compressive  load  (P)  and  end- 
shortening  (J}  Is; 


(P<P.) 


(2.21) 


(2.22) 


The  strain  energy  U of  the  entire  bean  is  the  sun  of 
strain  energy  Ui  in  the  cracked  portion  of  the  bean  and 
strain  energy  Vg  in  the  uncracked  portion  of  the  beam,  end 
it  has  the  following  form: 

£/  = £/,  +y, 

- f3/"o  4£/=s'']  P’(f,-a) 

Aa'  j 2E4 

(2  23) 

Strain  energy  release  rate  (6)  at  Che  crack  tip  is  obtained 
by  differentiating  total  strain  energy  tl  with  respect  Co 
half  of  crack  length  a [53] - Then  the  expression  for  G can 
be  derived  as: 


C = - 


2B 


I f ^ SE/!t'S(Sat2l) 

■2«|(4o  + 2i)’  (4a*2aLy 


. I6a'«*(12a*t4af,)  I2E/V1 
A(4a’ t2a'L)‘  Aa‘  | 


2.2.2  Delaminated  Seam  with  Throuoh-the-Thlckness  Stitches 


A delaminated  beam  with  stitching  is  depicted  in 


Fig-jre  2.5  A stitched  deleininated  heem 


First  let  us  consider  the  deleninated  portion  of  the 
with  stitching.  He  assume  that  the  buclcled  shape  of 
sublaminstes  is  simiiar  to  that  without  stitches 
described  in  the  previous  section.  Thus  we  assume  Us*cx 
w*tJsin^ffii/2a  where  unknown  coefficients  c and  d will 
different  for  the  unstitched  case.  The  deformed  shape 
the  cracked  portion  is  depicted  in  Figure  2.6, 


\ \ \ 


K n 

Figure  2.6  Deformed  shape  of  crack  portion  with  stitch 


Strain  energy  V of  the  delamated  portion  of  the  beam 
with  stitching  includes  two  parts:  strain  energy  Cliotucoh 
from  the  delaminated  portion  of  the  beam  without  stitching 
and  the  energy  contributed  by  the  stitches.  The  total 

strain  energy  0 is  of  the  following  form; 

tl  = U ....  . 

-V *]lK{2w)’dx  (2.25) 

here  K is  foundation  constant  Chat  represents  the  stiffness 
of  the  smeared  stitches.  The  foundation  constant  R is 
related  to  Che  actual  stitch  parameters  as: 


(2.26) 
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here  A.  is  the  area  of  cross  section  of  a stitch  (2  bobbin 
yarns;;  Em  is  Voung's  nodulus  of  stitch  aaterial  (bobbin 
yarn) t 2h  is  the  total  of  thickness  of  uncracked  beam;  N is 
the  stitch  density  (number  of  stitches/unit  area),  and  B is 
the  width  of  the  beam  in  the  y-direction. 

After  using  a procedure  similar  to  that  for  an 
unstitched  beam,  we  derive  the  buckling  load  P^t  of  a 
delaminated  beam  with  stitches  as  follows: 

p - ^ 3Jr(2g)  27) 

'■(2o)' 

From  Equation  2.2?  we  can  see  that  the  first  part  of  Per  is 
the  same  as  the  buckling  load  of  an  unstitched  beam,  and 
the  second  part  of  the  equation  is  contributed  by  the 
stitches.  Therefore  stitching  indeed  can  improve  the 
buckling  load. 

Now  we  can  extend  the  result  to  the  entire  stitched 
beam.  The  deformation  of  the  beam  undergoing  post-buckling 


The  relation  between  conpreesive  loads  F and  end 


shortening  6 Is: 
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When  compressive  loads  P Is  greater  than  the  buckling 
load  Pe,  of  the  stitched  beam  total  strain  energy  V of  the 
entire  delaainated  beam  with  stitching  is; 


j,_P\L-a)  ^ 3/”o  a ^ l2Ka'  , 

2JiA  2EA  Ea''  a-  ' 


(2  30) 


strain  energy  release 


g=--L^ 

2fl  ai|. 


J..  ] fs76A:'(20o*  + 12ii'I)  ^ 192a(4o= 

'”“'2fl|  (4a  + 2i)’il4ff'  .4(4a  + 2i)= 


*»KSl»o'  * 6a’ 1)  480  720XV1 
»’(4a  + 2/.)’  ' A ' EAx‘  ) 


Let  US  define  the  nondimensional  parameters  as; 

Pc,'-Pc, /SA,  ff-S/Bh,  a'=a/L,  B'-B/L,  h'-h/L,6  '-S/L,  K'‘K/E.  Per' 
is  non-diinensionei  buckling  load;  S'  Is  the  nondimensional 
energy  release  rate:  a'  is  the  nondimensional  half  crack 
length  of  the  structure;  S'  is  the  ratio  of  the  width  of  the 
structure  to  the  half  length  of  the  beam;  h*  is  the  ratio  of 
the  thickness  to  the  length  of  Che  structure;  S is  the 
nondimensional  end  shortening  of  the  beam  along  the  axial 
direction;  If  is  Che  nondimensional  stiffness  of  the 
stitches.  After  tedious  algebraic  manipulations  we  can 
derive  the  following  relations: 


(2  32) 


2S'  >rV(6a'’  + 2a)  288/C''(2Qg* 

■(4o  t2)’'^  9(4o'’  + 2a‘)’  «'B'A’(4a  *2)' 

-I-:)  8y*(4gU4o)  IKh 
* 3(4o“  + 2o')'  ' S'(4o'  + 2)’  ’ B 


24BB(8a’^-6o’)  ir'h'  360B’g' 
»’fl»(4<i  +2)“  '24a'*  a‘S'h‘ 


In  the  following  section  the  results  of  the  analytical 
model  are  discussed  in  term  of  the  nondimenslonal 
parameters . 

2.3  Results  and  Discussion 

The  results  presented  In  this  section  are  based  on 
Equation  2.32  and  2.33.  In  all  numerical  examples  h'-0.02 

Figure  2.6  shows  the  variation  of  buckling  load  Per’  as 
a function  of  semi-crack  length  a'  for  various  values  of 
stitch  parameter  K.  From  Figure  2.B  we  can  observe  that  the 
non-dimensional  buckling  loads  P,;e'  is  not  affected  by 
stitching  when  the  range  of  the  nondimensional  crack  length 
a' is  from  0.05  to  0.1.  When  a' is  greater  than  0.1,  P„'for 


35 

an  unstitched  beam  decreases  as  the  nondimensional  crack 
length  a'  increases.  This  is  understandable  because  a 
structure  buckles  more  easily  with  increasing  crack  length. 
When  we  look  at  stitched  beams  in  the  range  0.Ka'<0.4  (see 
Fig.  2.9}  the  buckling  load  is  always  higher  than  that  of 
unstitched  beam,  and  further  the  buckling  load  P„’  Increases 
as  the  stiffness  of  stitching  If  increases.  In  a stitched 
beam  an  interesting  phenomenon  occurs;  the  buckling  load  Pc,' 
decreases  as  the  crack  length  a'  increases  in  the  range 
0.Kak0.2,  however  tor  0.2<a'<0.4  Che  buckling  load  Pc,’  does 
not  decrease  due  to  Increasing  crack  length  a',  but  in 
contrast  P,,  increases  Co  some  extent  depending  on  the 
stiffness  of  stitching.  The  reason  is  Chat  Che  longer 
crack,  Che  more  stitches  involved  in  the  structure  and  the 
stitches  hold  the  structure  together  more  tightly, 
therefore  the  buckling  load  P,,’  can  be  improved  in  Che 
stitched  beam  even  with  longer  delamination.  A similar 
behavior  was  obtained  by  Shu  and  Mai  [39]  for  the  case  of 
thin  film  delamination  with  stitches.  Figure  2.9  shows  this 
phenomenon  in  detail. 

In  Figure  2.10  the  relation  between  nondimensional 
energy  release  rate  S' and  nondimensional  crack  length  a' and 
the  effect  of  K'  on  energy  release  rate  S'  for  a 


presented. 


nondimensionsl  end  shortening  5’»0.1  are 
crack  length  a'  is  less  than  0.1  the  stitching  has  no  effect 
on  non-dintenaional  energy  release  rate  G* . As  the 
nondimensional  delamination  size  a' grows  beyond  0.1  the 
effect  of  stitching  becomes  significant.  The  energy  release 
rate  s' decreases  in  a stitched  beam  compared  to  that  in  an 
unstitched  beam.  The  larger  the  stiffness  of  the  stitching, 
smaller  is  the  nondimensional  energy  release  rate  in 
stitched  beams.  Therefore  using  a stitching  can  decrease 
the  energy  release  rate  at  the  crack  tip,  and  the  crack 
propagation  can  be  delayed.  The  crack  will  propagate  when 
the  energy  release  rate  is  greater  than  the  critical  energy 
release  rate  (fracture  toughness  (Sc).  The  variations  of  G' 
have  similar  tendency  in  both  stitched  beam  and  unstitched 
bean.  Nondimensional  energy  release  rate  G'  dramatically 
increases  with  crack  Length  for  smaller  delaminations,  and 
reaches  maximum  G.„'  at  certain  crack  length  a'  and  then  O' 
goes  down  as  the  delamination  size  increases.  Whether  a 
crack  will  propagate  or  not  depends  on  fracture  toughness 
Gc.  Whether  it  will  be  a stable  crack  propagation  or  an 
unstatble  propagation  depends  on  the  rate  of  change  of 
energy  release  rate  G with  delamination  size  a.  Referring 
to  Figure  2,10,  we  can  note  that  in  the  beginning  the  crack 


37 

will  propagate  rapidly  in  an  unstable  manner.  Por  a’>0.1, 
dG'/cIa'  la  negative  and  hence  the  crack  propagation  will  be 

In  Figure  2.11  the  relation  between  nondimensional 
energy  release  rate  G'  and  nondimensional  stitch  stiffness 
at  a given  crack  length  |a'*0.2|  is  examined  and  the  effect 
of  non-dimensional  end  shortening  5'  is  shown.  It 
demonstrates  that  energy  release  rate  C'  decreases  with 
increasing  the  stiffness  K'  of  stitching.  The  more  non- 
dimensional  end  shortening  6'  is,  Che  more  nondimensional 
energy  release  rate  G'  has. 

The  variation  of  G'  with  S'  for  various  stitch 
stiffnesses  is  shown  in  Figure  2.12.  The  results  correspond 
to  a'*0.2.  The  relations  were  very  similar  for  other  crack 
lengths  also.  From  Figure  2.12  one  may  note  that  the  G' 
remains  close  Co  zero  in  the  beginning,  but  increases  after 
certain  critical  S'.  The  variation  of  G'  is  approximately 
quadratic  in  <5'.  Further,  it  may  be  noted  that  for  a given 
S',  G'  decreases  with  the  ft’.  This  is  brought  out  clearly  in 
Figure  2.12.  where  G'  is  plotted  as  a function  of  K'  for 
various  S'.  One  can  see  that  the  dramatic  decrease  in  G’ with 


increasing 


In  sumioary  stitching  always  increases  the 
load  for  buckling  in  delaminated  beams.  It  reduces  the 
energy  release  rate  at  the  crack-tip  thus  preventing  cr 
delaying  the  onset  of  delaminatlon  propagation.  Both  these 
effects  are  very  much  dependent  on  the  stitch  parameter  K' 
which  Is  given  by  (see  Equation  2.26) i 


K 


A,E,  NB 
2/t  E 


(2.34) 


From  Equation  2.34  one  can  see  that  K'  is  proportional  to 
the  stitch  axial  stiffness  and  stitching  densityf  but 
inversely  proportional  to  the  equivalent  Young's  modulus  of 
the  parent  laminate.  Thus  one  needs  to  use  high  density. 


high  stiffness  stitches  for  stiffer  laminates. 


Figure  2.8  Variation  of  with  a'  for 
in  the  range  0<a'<0.4 


different 


Figure  2.5  Vatiacion  of  with  a' 
for  different  K' 
in  the  range  0.1<a'<0.4 


Figure  2.11 


for  s given  a' 


43 


a given  a'  (a'-0.2) 


Figure  2.12  Variation  of  G’ 


S'  for 


CHAPTER  3 

FINITE  ELEMENT  ANALYSIS  OF  A DELAMINATED  CCMPOSITE  WITH 
STITCHING 

As  discussed  in  Chapter  2 delaoinatlons  can  greatly 
reduce  the  conpressive  load  carrying  capacity  of  laminated 
composite  structures.  Through-the-thickness  stitching  seems 
to  be  helpful  in  increasing  the  critical  load  at  which  a 
delaminated  beam  becomes  unstable  under  coi^ressive 
loading.  In  Chapter  2 an  analytical  solution  for  post- 
buckling  of  a stitched  composite  with  a delamination  was 


closed  form  solution.  In  this  chapter  finite  element 
analysis  is  used  to  investigate  the  effect  of  stitching  on 

Although  analytical  models  provide  insight  into  the  effects 
of  stitching,  more  detailed  behavior  of  a structure  can  be 
obtained  using  finite  element  models.  Several  restrictions 
of  Che  analytical  approach  can  be  easily  removed  in  FE 

is  homogeneous  and  orthotropic.  Whereas  FE  models  can 
handle  laminated  composites.  Unlike  in  the  analytical 
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Method  no  restriction  is  pieced  on  Che  location  of  Che 
craclc  in  the  thickness  direction.  Further  the  stitches  do 
not  have  to  be  smeared  as  continuous  springs;  each  stitch 
can  be  modeled  individually/  if  necessary.  In  this  chapter 
the  FE  model  is  used  to  understand  Che  effects  of  stitch 
parameters/  especially  stitch  density  and  stitch  yarn 
material  properties/  on  the  load-end  shortening  relation 
and  the  energy  release  rate  of  a delaminated  con^osite 
under  axial  con^ressicn. 


3.1  Specimen 

The  specimen  considered  for  finite  element  analysis  is 
made  up  of  4B  unidirectional  AS4  graphite/epoxy  plies.  The 
material  properties  for  each  lamina  are:  £ib19.4X10‘  psi/ 
Epl.3X10‘ psi/  Sir*0.1X10‘  psi.  i-i,-0.2a.  The  length  i of 
the  specimen  is  S inches  (127mm) . The  width  S is  1 inch 
125.4mm)  and  the  Chlc>:ness  H is  0.27  inch  (6.86mm)  . A 
single  delamination  (length  of  2a)  is  located  at  a distance 
h from  the  top  surface  of  the  beam.  The  specimen  is 
depicted  in  Figure  3.1.  The  specimen  is  under  axial 
compression.  The  FB  analysis  is  used  to  simulate  the 
compression  test.  In  particular  we  are  interested  in  the 
load-end  shortening  relationship  and  energy  release  rate  as 


function  of 


considered:  Glass 


Three  stitch  yarn  materials  were 
(3ST0  denier},  Kevlar  (16G0  denier]  and  Glass  2 (5952 
denier]  bobbin  yarn.  The  stitch  densities  are  4X1/4"  and 
0X1/8".  For  example  4X1/4"  stitching  pattern  means  a 
number  of  16  stitches  per  square  inoh  where  the  pitch  is 
1/4  inch  and  the  distance  between  two  adjacent  stitch  lines 


Figure  3.1  A delaminated  beam-like  structure  with  stitche 


3.2  Finite  Element  Model 

The  commercial  finite  element  software  ABAQUS  |54J  was 
used  to  simulate  the  stitched  composite  with  through  the 
width  delamination.  This  specimen  is  modeled  by  two 


beam  are  distributed  uniformiy  along  the  axial  direction. 
Type  B21  element,  which  conaiders  tranaverse  sheer  effect 
is  adopted.  The  stiffness  matrix  of  this  type  is  given  In 
ABAOUS  Theory  Manual  ISS].  The  element  size  is  chosen  as 
0.0125  inch.  In  an  undelaminated  portion  the  connection 
between  top  and  bottoms  is  aohieved  by  using  the  EQUATION 

depicted  in  Figure  3.2. 


The  relations  between  the  top  and  bottom  beam 


Jisplj 


4B 


(3  I) 

(3.2) 

(3.3) 


the  X direction,  w is  the  nodal  displacement  in  the  z 

Ihcough-the-thlckness  stitches  are  modeled  as  linear 
springs  using  SPRIN52  elements  in  ABAQUS.  Two  spring 

represents  the  resistance  offered  by  the  stitch  to  relative 
transverse  displacements  of  the  top  and  bottom  beams.  The 
second  spring  element  represents  the  relative  displacement 
in  Che  axial  (longitudinal)  direction.  These  spring 
elements  connect  the  two  corresponding  nodes  on  the  top  and 
bottom  beams.  The  relative  displacement  across  a SPRIHG2 
element  is  the  difference  between  the  1th  component  of 


ve  denote  the  axial  and  transverse  displacements  by  u and 
w,  then  relative  displacements  are  [see  Figure  3.3) : 

Ai;  = v, -V,;  d»-  = »,  (3  4) 


Figure  3.3  Node  specification 
for  SPBING2  element 


The  procedure  for  determining  the  spring  constants  of 
the  stitches  is  as  follows.  Figure  3.4  shows  a stitch  where 
the  beam  surfaces  have  relative  displacements  in  both  x and 
i directions.  The  initial  length  of  the  stitch  is  the  total 
thickness  of  the  beam,  H.  The  relative  displacements  in  the 
X and  a directions  are  du  and  dw.  respectively.  The 
elongation  of  the  stitch  dL  is  related  to  the  displacement 
components  by 


^ = &Lcos9 

Hi) 


several  buckling  mode  shapes.  Thus  a linear  bifurcation 
analysis  has  to  be  performed  to  determine  the  buckling  mode 
shapes.  In  Che  present  analysis  only  the  first  mode  shape 
was  considered  In  Che  Iniciel  Imperfection.  The 
imperfection  amplitude  was  51  of  the  eigenoiode  calculated 
by  ABAQUS. 


Top  beam 


Figure  3.4  Stitch  deformation 


The  nonlinear  analysis  used  the  Riks  algorithm  in 
order  Co  determine  the  successive  equilibrium  positions.  To 
analyte  a post-buckling  problem,  it  must  be  turned  into  a 
problem  with  continuous  response  instead  of  bifurcation. 
This  effect  can  be  accomplished  by  introducing  an  initial 
imperfection  into  a "perfect"  geometry  so  that  there  is 
some  response  in  the  buckling  mode  before  Che  critical  load 
is  reached.  The  post -buckling  problem  is  a geometrically 
nonlinear  static  problem.  Static  egullibrium  states  during 


the  unstable  phase  of  the  response  can  be  found  by  using 
Che  modified  Biks  method.  This  method  is  used  for  the  case 
where  the  load  magnitudes  are  governed  by  a single  scale 
parameter.  The  method  can  provide  solutions  even  In  the 
case  of  oomplex/  unstable  response.  The  Biks  method  uses 
the  load  magnitude  as  an  additional  unknown;  it  solves 
simultaneously  for  load  and  displacements.  Therefore, 
another  quantity  must  be  used  to  measure  the  progress  of 
Che  solution;  ABAQUS  uses  the  arc  length,  1,  along  Che 
static  equilibrium  path  in  load-displacement  space.  This 
approach  provides  solutions  regardless  of  whether  the 
response  is  stable  or  unstable.  The  modified  Riks  algorithm 
is  described  in  detail  in  [55] . 

Energy  Release  Rate 

The  purpose  of  stitching  or  any  other  forms  of 
through-the-thickness  reinforcement  is  to  increase  the 
buckling  load  and  also  Co  prevent  Che  cracks  from 
propagating  further.  Energy  release  rate  at  the  crack  tip 
has  been  found  to  be  a useful  parameter  in  predicting  crack 
propagation  in  laminated  composites.  The  delaminacion  is 
expected  to  propagate  further  when  the  energy  release  rate 
exceeds  the  corresponding  fracture  toughness  value  of  the 
composite  material  system.  DelaminaCions  have  a higher 


propensity  for  propagstion  along  the  plane  of  delamlnation 
as  the  interlaminar  fracture  toughness  is  much  lower  for 
most  laminated  composites.  In  particular  Mode  1 (opening 
mode)  interlaminar  fracture  toughness  is  about  one-halt  of 
Mode  II  (shearing  model  fracture  toughness.  For  example, 
the  values  for  a typical  AS4/3501-6  graphite/epoxy 
composite  are  [46] : Gje  * 300  N/m  and  * 6*^0  N/m.  Thus 
the  effectiveness  of  stitching  can  be  judged  by  the 
reduction  in  energy  release  rate,  G,  at  the  crack  tip. 

There  are  three  different  methods  to  compute  G:  (il 
computing  stress  intensity  factor  at  the  crack-tip;  (lil 
Strain  Energy  Density  Method;  (lii)  Virtual  crack  growth 
method.  The  first  method  requires  detailed  three- 
dimensional  stress  field  near  the  crack  tip,  which  can  only 
be  obtained  by  using  three-dimensional  finite  elements  with 
a fine  mesh.  Since  we  ace  using  structural  elements  in  the 
present  analysis,  the  stress  intensity  factor  approach 
cannot  be  used.  The  strain  energy  density  method  developed 
by  Sankar  [521  requires  computation  of  energy  densities 
just  behind  and  ahead  of  the  crack  tip.  However  in  the  case 
of  stitched  composites,  the  stitch  forces  affect  the  energy 
densities  significantly,  and  hence  a very  fine  mesh  is 
needed  near  the  crack  tip-  Further  ABACUS  does  not  provide 
the  energy  density  values  readily,  and  they  need  to  be 
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computed  from  the  force  and  moment  resultant  Just  ahead  and 
behind  the  cracfc  tip.  Hence  the  third  method  is  used  In  the 
present  study,  and  it  is  explained  below. 

Energy  released  rate  G is  computed  by  using  the 
following  formula  153): 


C = - 


2B  da 


(3  8) 


here  G is  the  strain  energy  in  the  beam  for  a given 
displacement,  a is  the  crack  length,  and  B is  the  width  of 
the  beam.  In  this  method  the  load-end  shortening  relation 
P-S  is  developed  using  the  aforementioned  nonlinear  FE 
analysis.  For  a given  displacement  (end  shortening)  the 
area  under  the  P-S  diagram  up  to  the  displacement  represents 
the  strain  energy  stored  in  the  beam.  Let  it  be  equal  to 
UiSi . the  change  in  strain  energy  is  computed  numerically  by 
repeating  the  nonlinear  FE  analysis  for  a crack  length 
at-Aa,  where  Aa  is  the  incremental  crack  length.  The  area 
under  the  new  P-S  curve  represents  the  strain  energy  for  the 
new  crack  length.  Then  G is  obtained  by  numerically 
differentiating  UiS)  (see  Fig.  3.5): 
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Figure  3.5  Strain  energy  change 
due  to  crack  growth 


3.5  Results  and  Discussion 

In  this  section  the  results  of  the  nonlinear  FE 
analysis  are  presented.  The  effects  of  stitching  are 
understood  by  analyzing  the  load-end  shortening  dlagrans, 

quantity  of  interest  is  the  energy  release  rate  S at  Che 
crack  tip.  The  S is  also  plotted  as  a function  of  S tor 
various  stitches.  Another  Important  and  interesting 
information  is  Che  amount  of  axial  strain  in  the  stitches, 


:ially  during 


:kling  process.  Although 


stitches  prevent  delaminatlon  propagation  by  reducing  the 
G,  they  ere  subjected  to  strain  by  the  beam  segments  that 
tend  to  buckle.  Thus  the  strains  in  the  stitches  are  also 
plotted  as  a function  of  S. 

example;  Kevlar  (1600  denier),  Glass  1 OSIO  denier),  and 
Glass  2 (59S2  denier).  Their  properties  are  listed  in  Table 
3.1.  Two  types  of  stitch  patterns  were  considered:  4x1/4" 
(16  stitches  per  square  Inch)  and  8xl/8"(64  stitches  per 

specimen  dimensions  were  described  in  section  3.1.  The 
delaminatlon  was  assumed  to  be  at  a distance  equal  to  l/4th 


The  P~S  dla^rame  for  16  ssi  (stitches  per  square  inch) 
and  6A  ssi  are  given  in  Figure  3.6  and  3.7,  respectively. 

It  is  clearly  seen  that  stitching  can  increase  the  peak 
load  by  20%'30%  depending  on  Che  stitch  eiaCerial  and  stitch 
density.  It  is  also  seen  that  the  difference  in  ataximuiQ 
loads  due  to  different  stitch  materials  is  not  significant. 
The  effects  of  stitch  density  on  Che  maximum  load  can  be 
inferred  from  Figure  3.6.  It  may  be  seen  that  there  is  a 
slight  increase  in  the  maximum  load  due  to  an  increase  in 
stitch  density  from  16  ssi  to  64  ssi.  The  results  presented 
in  Figure  3.6  through  3.8  are  for  a crack  length  of  1 inch. 

Figure  3.9  shows  the  effects  of  various  stitches  on 
the  maximum  load  as  a function  of  the  crack  length,  It 
should  be  noted  that  the  beam  length  is  5 inches,  however 
the  effective  length  is  only  3 Inches  (see  Figure  3.2). 

From  Figure  3.9  it  is  seen  that  stitching  has  no  effect  for 
very  short  delamlnatlons  (a<l  in).  For  O.S  in.  < a < 1.5 
in.  the  64  ssi  specimens  carry  almost  the  same  maximum  load 
as  the  undelaminated  specimens.  That  is,  the  high  density 
stitching  is  able  to  completely  suppress  delamination 
buckling  and  restore  the  original  load  dairying  capacity. 
The  16  ssi  specimens  also  show  marked  improvement  in  the 
maximum  load.  However  for  a>2.0  inches  the  effect  of 
stitching  seems  to  diminish.  This  is  because  the 


delaminacion  length  approaches  the  total  length  of  the 
beam,  and  global  buckling  takes  over  the  sub-delamination 
buckling.  During  global  buckling  the  delamination  is  under 
more  shearing  mode  than  opening  mode,  and  Che  stitches  are 
not  that  effective  in  suppressing  the  relative  sliding 
notion  of  the  delaminated  surfaces. 

The  effects  of  stitching  on  energy  release  rate  G are 
shown  in  Figure  3.10  through  3.12.  Figure  3.10  corresponds 
to  a stitch  density  of  16  ssi.  It  is  seen  that  G is  almost 
negligible  until  the  sublamlnate  buckling  begins. 

Thereafter  the  6 rises  rapidly  with  the  end  shortening  of 
the  beam.  The  stitches  are  able  to  delay  the  point  when  G 
begins  to  rise.  For  example,  in  Figure  3.10,  G begins  to 
rise  at  5«0.03  in.  in  an  unstitched  beam,  whereas  it  is 
delayed  until  6«0.55  in.  in  stitched  beams.  Further,  the  G 
is  greatly  reduced  compared  to  the  unstitched  beam.  For 
5»0.1  in.  Che  value  of  G is  reduced  from  about  600  lb. /In. 
in  an  unstitched  beam  to  about  100  - 200  lb. /in.  in 
stitched  beams.  Further  it  can  be  seen  that  stitches  with 
higher  axial  rigidity  (see  Table  3.1)  cause  the  greatest 
reduction  in  G. 

Similar  results  for  64  ssi  stitching  are  presented  in 
Figure  3.11.  Both  stitch  densities  are  compared  with 


unstitched  beam  in  Figure  3.12.  It  nay  be  seen  that  64  ssl 
stitching  offers  slightly  higher  reduction  in  S. 

As  mentioned  earlier  monitoring  the  strains  in  the 
stitches  is  an  important  consideration  in  analyzing 
stitched  specimens.  Because  stitch  yarns  could  break  if  the 
strains  exceed  the  allowable  limit  and  this  could  trigger  a 
domino  effect  of  successive  stitch  failures  and  dynamic 
delamination  propagation^  and  lead  to  catastrophic  failure. 
The  strains  in  various  stitches  as  a function  of  end 
shortening  are  shown  in  Figure  3.13.  It  should  be  noted 
that  the  crack  length  is  1 inch  and  the  stitch  pattern  is 
4x1/4".  Stitch  3 (Figure  3.13)  is  at  the  center  of  the  beam 
and  thus  undergoes  large  strains.  Stitch  1 is  actually  at 
the  crack  tip  and  hence  there  is  no  strain  at  all.  Stitch  2 
is  in  between  stitches  1 and  3.  The  variation  of  stitch 
strain  with  end  ahortening  is  strikingly  similar  to  that  of 
G.  discussed  earlier. 

The  strains  of  stitches  located  at  the  center  of  beam 
in  various  stitch  materials  and  stitch  patters  are  shown  in 
Figure  3.14.  It  may  be  seen  that  the  strains  are  sensitive 
to  the  stitch  material  as  well  as  the  stitch  density.  In 
general  the  strains  ate  significantly  less  for  higher 
stitch  density.  Further,  stitches  with  lower  axial 
stiffness  (A.S.)  undergo  larger  strains. 
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From  the  above  discussion  we  can  arrive  at  the 
following  conclusions: 

1.  Stitching  can  improve  the  buckling  loads  of 
delaminated  beams.  The  increase  is  not  sensitive 
to  the  three  types  of  yarns  used  in  this  study. 

The  yams  had  axial  rigidities  in  the  range  of 
13,000  S - 223,000  N and  they  seem  to  perform 
equally  well.  However,  Che  maxiinuin  load  increases 
with  Che  stitch  density.  This  is  consistent  with 
the  experimental  observation  of  Sharma  and  Sankar 
146] . 

2.  Through-the-thickness  stitches  significantly 
reduce  the  energy  release  rate  at  Che  crack  tip. 
They  also  increase  Che  end  shortening  required  to 
cause  sublaminaCe  buckling  when  the  energy  release 
rate  rises  rapidly.  By  reducing  Che  energy  release 
rate  stitching  prevents  Che  crack  propagation 
which  may  lead  Co  catastrophic  failure  of  the 
specimen.  The  axial  rigidity  of  the  stitches  also 
affects  Che  energy  release  rate.  The  Higher  axial 
rigidity  the  greater  is  the  reduction  in  G. 

3.  The  maximum  strain  in  the  stitches  depends 
strongly  on  Che  axial  rigidity  of  the  stitches  and 
Che  stitch  pattern.  The  stitches  with  small  axial 


rigidity  undergo  large  axial  straina  during 
sublaninate  buckling. 


Figure  3-6  Compressive  load  vs.  end  shorten 
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Figura  3-7  Compressive  load  vs.  end  shortening 
at  crack  length  (a*I  in,) 


Figure  3.8  Compressive  load 


shortening 


Crack  length  a (in. ) 


End  shortening  5 (in.) 

Figure  3.10  Energy  release  rate  vs.  end  shortening 
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End  shortening  2 (in.) 
figure  3.11  Energy  release  rate 
vs. end  shortening  (a*l"  h/H^l/4) 


End  shortening  5 (in.) 


figure  3.12  Energy  release 


shortening 


End  shocCsnlng  6 (in,) 


Figure  3.13  Stitch  etraine  vs.  end  shortening 
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End  shortening  6 fin.) 


Figure  3.14  Maximum  strain  in  a stitch 
vs.  end  shortening  ta-l"  h/H-1/4) 


stitch  ir 
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The  impact  problem  is  depicted  in  figure  4.1.  The  force 
acting  on  the  beam  is  not  known  a priori;  and  has  Co  be 
calculated  as  part  of  the  solution  of  the  problem.  A major 
difference  between  Che  impact  loading  and  the  compressive 
loading  discussed  in  the  earlier  chapters  is  Chat  the 
delamination  is  under  pure  Mode  II  (shearing  mode) 
conditions.  Previous  experimental  studies  b/  Sharma  and 
Sankar  [46]  have  found  that  under  pure  Mode  II  Che  stitches 
try  to  plough  through  the  matrix,  and  the  resistance  offered 
by  the  matrix  is  responsible  for  Che  Increase  in  apparent 
fracture  toughness.  Thus  Che  stitch  model  has  to  be  modified 
to  account  for  this  phenomenon. 


Impactor 


* 


M,:  Impactor  mass 
V,!  Impact  velocity 


t:  Time  Variable 


Figure  4.1  A structure  under  impact 


In  the  following  sections  we  discuss  Che  assumptions 
made  in  simplifying  the  impact  problem,  Che  stitch  model  and 
the  impact  simulation  procedure.  The  results  focus  on  the 
extent  of  delamination  propagation  due  to  impact  in 
specimens  with  different  stitch  densities. 

4.1  Basic  Assumptions 

The  following  assumptions  are  made  to  simplify  the 
impact  simulation: 

1.  The  velocity  of  impact  is  low  compared  to  the 
velocities  of  wave  propagation  in  the  composite 

2.  The  projectile  is  assumed  to  be  rigid  compared  to 
the  target.  Therefore  the  impactor  can  be  treated  as 
a rigid  body  and  its  equation  of  motion  is  greatly 
simplified. 

3.  The  target,  laminated  beam  in  the  present  case,  is 
highly  flexible.  The  deflection  of  the  beam  is 
expected  to  be  much  higher  than  the  local 
Indentation,  and  hence  the  Hertzian  indentation 
effects  can  be  neglected. 

4.  The  impactor  mass  is  much  greater  that  that  of  the 
beam,  and  hence  the  impact  duration  will  be  very 
long  compared  to  the  fundamental  period  of  vibration 
of  the  beam.  Therefore  the  target  can  be  represented 


by  a simple  spring-mass  system.  The  stiffness  of  the 
spring  can  be  approximated  by  the  static  beam 
stiffness  k at  the  impact  location. 

5.  An  embedded  central  delamination  is  assumed  to  be 
pre-existing  in  the  laminated  beam.  The  delamination 
is  symmetrically  located  in  the  simply  supported 

Because  of  Assumption  S above  we  need  to  compute  only 
the  static  stiffness  of  the  delaminated  beam  due  to  a 
central  transverse  force.  We  will  assume  that  the 
delamination  will  propagate  along  the  same  plane  when  the  G 
exceeds  Gr/i,,  Mode  II  fracture  toughness  of  the  parent 
laminate . 

The  resistance  offered  by  the  matrix  under  shearing 
mode  is  computed  as  follows.  The  ploughing  resistance  of  the 
stitches  can  be  represented  as  a distributed  shear  traction 
(force/unit  length)  at  the  interface  of  the  two  sublaminates 
that  are  stitched  together.  The  traction  pn  is  estimated  as: 


here  B is  the  beam  width,  D is  the  diameter  of  stitch  yarn, 
h is  the  greater  of  the  two  sublamlnate  thic)cnesses,  Oy  is 
the  yield  stress  of  the  surrounding  matrix,  s Is  the  stitch 
spacing  in  the  width  direction,  and  n is  the  number  of 


stitches  per  inch.  In  terms  of  stitch  density  the 
distributed  traction  can  be  written  as; 


A = NBDha, 

where  N*l/(nKS)  is  the  number  of  stitches  per  unit  area.  As 
the  delamlnatlon  propagates  new  stitches  come  into  action  in 
the  freshly  created  delamlnatlon  areas,  and  they  offer 
additional  shear  resistance.  This  assumption  is  consistent 
with  the  experimental  observations  of  Sharma  and  Sankar 
ME). 


e.2  Analytical  Model 

The  equation  of  motion  of  the  impactor  along  with  the 
initial  conditions  can  be  written  as: 


(4. 1 -a) 
(4.1-i) 


here  Mo  is  the  impactor  mass,  q Is  the  impactor  displacement 
which  is  same  as  the  transverse  deflection  of  the  target 
beam  at  the  point  of  impact,  is  the  initial  velocity  of 
the  impactor  or  the  impact  velocity,  and  F(ql  is  the  contact 


force.  The  contact  force  F is  a function  of  the  beain 
deflection.  The  contact  force  F will  be  a linear  function  of 
q if  there  were  no  stitches  or  deiamination  propagation. 
However  in  the  present  case  it  will  be  a nonlinear  function. 
Once  F(q}  is  determined,  then  the  equation  of  motion 
Equation  4.1a-4.1-c  can  be  numerically  integrated  to  obtain 
qttl.  From  q(ti  one  can  compute  the  impact  force  history 
Fit)  using  the  F-q  relation.  In  the  following  section  we 
discuss  the  procedure  for  determining  the  F-q  relation  in  a 
stitched  delaminated  beam. 

^.2.1  Relation  between  Contact  Force  and  Beam  Deflection 

The  problem  to  be  solved  in  this  section  is  depicted  in 
Figure  4.2.  Because  of  symmetry  only  one  half  of  the  beam 
will  be  analyzed.  The  problem  is  to  find  the  relation 
between  transverse  force  F and  deflection  q at  the  center  of 
the  beam  in  the  delaminated  stitched  bean.  Further  the 
energy  release  rate  G at  the  craclc  tip  needs  to  be  computed 
also.  In  the  numerical  simulation  the  cracb  will  be 
propagated  by  a small  distance  (symmetrically  on  both  sides) 
if  the  5 exceeds  the  Mode  II  fracture  toughness  G,ic  of  the 
parent  laminated  material  system. 

First  we  will  provide  an  overview  of  the  procedures  to 
be  followed.  Since  the  structure  is  symmetric  we  can  analyze 


L Is  imcrac)ced  length  in  Figure 


length  is  denoted  by  e and 
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Figure  4.2  K delaminated  bean  with  stitches 
subject  to  impact  load  F 


Figure  4.3  Right  half  side  of 
specimen  in  Figure  4.2 


The  right  half  beam  can  be  divided  into  three  elements: 
Element  1,  Element  2 and  Element  3.  In  Figure  4.3  they  ace 


body  dlsgi 
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denoted  as  <1>,  <2>  and  <3>.  Free 
three  elements  are  shown  In  Figure  A,i.  At  each  end  of  the 
element  there  are  three  displacements,  u,  u,  and  v-  and 


represents  the  couple  corresponding  to  the  rotation  i/>.  We 
use  the  Timoshenko  beam  theory  in  modeling  the  beams.  The 
equations  of  equilibrium  of  a general  element  are  given  in 
the  l^spendix.  These  equations  can  ha  easily  solved  as  shown 
in  the  Appendix  to  obtain  a relation  between  the  six  forces 
acting  at  the  two  ends  of  the  element  and  the  six 
corresponding  displacements.  In  order  to  relate  the  forces 
in  the  three  different  elements  we  use  force  and  moment 
equilibrium  equations,  and  compatibility  equations  at  the 
junctions  (nodes)  where  the  elements  meet.  Further,  the 

implemented.  This  procedure  is  similar  to  assembling 
elements  in  the  Finite  Element  Method  except  the  stiffness 
matrix  of  each  element  is  obtained  exactly  by  solving  the 
differential  equations  of  equilibrium.  The  forces  and 


displacements  in  each  element 


displacements  at  the  junction  of  all  three  elements  is 
depicted  in  Figure  4.3. 


Figure  4.4  Free  body  diagrams  of  part  <1>,  <2>,  <3> 


Figure  4.5  Compatibility  of  displacement 
in  axial  direction  at  joint  position 


The  force  and  moment  equilibrium  equations,  Che 
coB^atibllity  equation  and  boundary  conditions  are  as 


follows: 


Force  end  moment  equilibrium  equetlone: 
F„  +F^  +F^  =0 
F„  +f„  =-- 
F„  *F,^  *F„  =0 

c.+c.+r,*/-- 

’ "2  “2 

Compatibility  equations  at  the  joints: 


(4  2) 

(4.3) 

(4.4) 

(4.5) 


(4.«) 


(4.7) 

(4.8) 


Boundary  conditions: 


(4.9) 


B,=«,=0  (4,10) 

= 0 (4.11) 

(4.12) 

“'»  = = C,  = 0 (4.13) 


ely.  Cj  is 


Odd  Indices  i fi-1,3,5)  denote  the  left  end  node  of  each 
element  and  even  indices  U=2,4,61  correspond  to  the  right 
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using  the  procedure  explained  in  the  next  paragraph.  If  the 
value  of  G exceeds  the  Mode  II  fracture  toughness  Gi:a*  the 
length  of  the  delanination  is  Increased  by  a small  amount. 
The  numerical  value  of  the  extension  is  arrived  by  a trial 
and  error  method  until  G equals  Gus  for  Chat  deflection 
increment. 

In  order  to  conqsute  the  G we  used  the  strain  energy 
density  method  derived  by  Sankar  (S2J.  This  method  is  very 
much  suitable  for  the  present  model  as  the  force  and  moment 
resultant  ahead  and  behind  the  crack  tip  can  be  obtained  in 
closed-form  from  the  solution  of  the  differential  equations 
of  equilibrium  and  the  energy  density  values  thus  calculated 
does  not  have  any  numerical  errors.  Consider  the  three 
elements  surrounding  the  crack  tip  as  shown  in  Figure  4.3. 
There  ate  two  (Elements  1 and  2)  behind  the  crack  tip  and 
one  ahead  of  the  crack  tip  (Element  3).  The  G is  derived  asr 

‘ <4.15) 

here  l/i  represents  the  strain  energy  density  (strain  energy 
per  unit  length  of  the  beam),  the  superscripts  (1)  and  (2) 
denote  the  cross  sections  immediately  behind  the  crack  tip 
and  (3)  denotes  the  cross  section  immediately  ahead  of  the 
crack  tip  and  B 1s  the  width  of  the  beam  in  the  y-directlon. 


The  strain  energy  density  in  teroia  of  force  and  moment 
resultants  is  given  by: 


where  P,  H and  V are  the  axial  force,  bending  moment  and 
shear  force  resultants;  EA,  El  and  5A  are  the  equivalent 
axial,  flexural  and  shear  rigidities  of  Che  beam  cross 

^.2.2  Impact  Response 

After  computing  the  f-q  relation  for  a beam,  Che  impact 
equations  (Equation  4.1-cl  can  be  solved  numerically.  The  F- 
q relations  were  stored  in  a spread  sheet  program  (Excel™)  . 
The  expression  for  F(q)  on  the  RHS  of  Equation  4.1-a  can  be 
approximated  by  a linear  interpolation  in  each  small 
integration  step  as: 


where  q,  and  qi  are  Irapactor  initial  and  final  displacements 
and  f,  and  f-  are  initial  and  final  contact  force  in  each 


(4,16) 


(4.J7) 


increment. 


equation  4.1-a  the  relation  between  contact  force 


F and  velocity 


can  be  derived  as  follows: 


{4.18-0) 


where 


Substituting  Equation  4.1B-b  into  Equation  4.18-a  and  a new 
equation  Equation  4.18-c  relating  the  contact  force  Ffg;  and 
the  Impact  velocity  V can  be  obtained: 


(4.18-c) 


Substituting  Equation  4.1S-b  Into  Equation  4.18-c  and 
integrating  both  sides  of  the  equation,  the  relation  between 
V and  F(q)  can  be  derived  as: 


-!•,=  ) = -I /-(gVi, 


(4.18-rf) 


Substituting  Equation  4.17  into  Equation  4.18-d  the  velocity 
can  be  expressed  in  following  form 


= •jy'’ - - iKF,  * F.) 


(4.19) 


initial  and  final  velocities  in  each 


The  corresponding  impact  time  can  be  derived  using 
Equation  4.16-b  by  assuming  that  the  impact  velocity  in  each 
step  varies  linearly; 


here  qj  and  qt  are  in^actor  initial  and  final  displacements 
in  each  step.  Integrating  Equation  4.18-b  and  Chen 
substituting  Equation  4.20-a  we  derive  Che  following 
integral  equation: 


Integrating  Equation  4.20-b  an  expression  for  time  ts  can  be 
obtained  as: 


here  ti  and  t;  are  initial  and  final  time  in  each  step. 
Equation  4.20-c  provides  the  q-t  relation  for  Che  impact 
problem.  From  that  using  F-q  relations,  we  can  obtain  the  F- 
C relation  or  the  impact  force  history.  Since  we  know  the 
delamination  length  at  each  displacement  in  the  static 


(4.20 -0) 


(9-  -9i) 


(4.20-4) 


I =/  I (9; 


(4  20-c) 


problem,  i.e.,  e-q  relation,  we  can  translate  that  into  a-t 
relation,  and  thus  the  propagation  of  delamination  can  be 
followed. 

4,3  numerical  Example,  Results,  and  Discussion 

The  specimen  dimension  and  material  properties  used  in 
Che  impact  simulation  ace  as  follows:  initial  uncraclced 
length  L^li.2  mm,  initial  crack  length  afl*27.4  mm,  beam 
width  fi-25.4  mm,  equivalent  Young!  5 modulus  £^90.375  GPa, 
equivalent  shear  modulus  69^=6. 6GPa.  Stitching  material  is 
3570  denier  glass  yam.  Nine  different  examples  were 
studied.  In  these  examples,  the  impactor  mass,  impact 
velocity,  position  of  the  delamination  in  Che  thickness 
direction  and  the  Mode  II  fracture  toughness  G„cwere 
varied.  The  parameters  used  in  these  examples  are  listed  in 

In  Table  4.1  H Is  the  total  thickness  of  the  beam 
and  h,  is  the  distance  of  delamination  from  the  top  surface 
(impact  surface).  For  each  example  three  different  cases  - 
beam  without  stitching,  16  ssl  stitches,  and  64  ssi  stitches 
- were  considered.  Thus  a total  of  27  impact  slmulacions 
were  performed.  The  simulations  were  stopped  when  the 
contact  force  becomes  equal  to  sero  denoting  the  contact 
between  the  impactor  and  the  beam  has  ceased.  In  two  cases 
(Example  2,  no  stitches  and  16  ssi  stitches)  the  simulation 


delanination  propagated 


has  to  be  stopped  when  the 


energy) . 

From  the  static  analysis  we  obtain  the  load-deflection 
relation  (F-q)  and  the  delamination-deflection  relation  (a- 
ql . Sample  F-q  relations  for  Example  1 are  shown  in  Figure 
4.7.  The  F-q  relation  is  linear  until  the  crack  begins  to 
propagate.  It  may  be  noted  that  the  load  at  which  the  ctaclc 
begin  to  propagate  is  almost  the  same  for  all  three  cases 
(no  stitches,  16  ssi.  and  64  ssi  stitches).  After  that  the 
curves  ta)te  different  shapes  depending  on  the  stitch 
density.  The  maximum  load  that  the  beam  can  carry  very  much 
depends  on  the  stitch  density.  The  64  ssi  beam  carries  about 
SOI  more  load  than  the  unstitched  beam.  The  unloading  was 


unloading  curve 


straight  line  joining  the  point  of  unloading  and  the  origin. 
This  assun^tion  is  validated  by  the  Mode  tl  experiments 
conducted  by  Shame  and  Sankar  [46]. 

Another  Interesting  result  that  can  be  deduced  from  the 
static  load-deflection  curve  is  the  apparent  fracture 
toughness  of  stitched  laminates.  The  area  enclosed  by  the 
load-deflection  diagram  (Fig. 4. 7)  denotes  the  work  of 
fracture.  Since  we  know  the  extent  of  delaminetion 
propagation  we  can  compute  the  apparent  fracture  toughness 


here  AH  is  the  work  done  and  AA  is  the  new  delamination 
surface  created.  The  apparent  fracture  toughness  for  various 
cases  is  presented  in  Table  4.2  along  with  that  tor 
unstitched  laminates.  The  numbers  in  parentheses  are  the 
percentage  increase  In  apparent  fracture  toughness.  It  may 
be  seen  that  the  percentage  increase  in  Gne  is  higher  for 
laminates  with  lower  fracture  toughness. 

The  results  for  each  Impact  analysis  include  the 
complete  impact  force  history  (F-tl  and  the  delamination 
propagation  history  la-t).  A sample  Impact  force  history  is 
shown  in  Figure  4.8.  In  general,  stitched  beams  carry  more 
impact  force,  provided  the  impact  energy  is  sufficient  to 
cause  delamination  propagation.  For  low  iiigsact  energies,  the 


impact  force  history  will  be  identical  in  stitched  and 
unstitched  beams,  because  the  stitches  come  into  effect  only 
when  there  is  sufficient  energy  to  propagate  Che 
delamination. 


£ab^^^2Com£arisons_of_Sii, 


Giic  for  16ssi 


The  reeulte  presented  in  Table  4.3  show  for  each  case 
impact  energy,  Gu,  of  the  parent  material  system,  the 
contact  force  F.  at  which  the  delamination  began  to 
propagate,  Che  final  crack  length  and  Che  maximum 

contact  force  f„.  during  the  impact  event.  The  extent  of 
delamination  propagation  is  also  shown  In  the  bar  charts  in 
Figures  4.9  and  4.10.  Figure  4.9  considers  the  examples 
wherein  the  G,/c  * S30  J/m*'  and  Figure  4.10  corresponds  to 
G;„=300  J/m*. 

There  ere  many  interesting  observations  that  can  be 
made  from  the  results  presented  in  Table  4.3.  The  contact 
force  at  which  the  delamination  propagates  is  almost  the 
same  in  unstitched  and  stitched  beams.  This  initiation  force 
depends  only  on  the  6,,c  and  the  position  of  the  crack 


(hl/H).  Thus  it  does  not  depend  on  impact  parameters  such  as 
impact  energy. 

In  general  the  extent  of  crack  propagation  at  the  end 
of  the  impact  event  is  the  least  in  the  64  ssi  beam  and  the 
highest  in  the  unstitched  beams.  The  result  for  16  ssi  beams 
are  somewhere  in  between.  This  can  also  be  observed  readily 
from  the  bar  charts  in  Figures  4.9  and  4.10.  However  the 
amount  of  delamination  propagation  depends  also  on  the 
impact  energy,  Giu  and  hl/H.  Stitching  is  very  effective  in 
the  beam  with  lower  inherent  fracture  toughness.  In 
Example  2 (Gns  * 300  N/m)  the  crack  propagates  all  the  way 
to  the  ends  of  the  beam  in  the  unstitched  and  16  ssi 
specimens,  whereas  the  crack  propagated  up  to  a=S9.4  mm  in 
the  64  ssi  beam.  In  Example  1 (6,«  = 530  N/ml  the  stitches 
were  able  to  reduce  the  delamination  extension  by  about  6 - 

Stitches  are  also  more  effective  when  the  crack  is  in 
the  middle  plane  of  the  laminate  (Al/N-O.S)  compared  to  the 
oases  wherein  the  crack  is  near  the  top  surface  of  the  beam 
.25) . This  can  be  explained  as  fallows.  In  an 
undelaminated  beam  the  shear  stresses  are  higher  at  the 
midplane  (1.5  times  the  average  shear  stress)  compared  to 
the  plane  at  1/4  distance  from  the  top.  Thus  when  the 
delamination  Is  at  the  midplane  the  tendency  for  propagation 
is  much  higher.  In  fact  the  energy  release  rate  G is  higher 


tor  midplane  delaminations-  Thus  the  stitches  play  a very 
useful  role  in  preventing  crack  propagation.  This  situation 
is  similar  to  Che  effectiveness  of  stitches  for  various 
G„c'S. 

The  stitches  become  more  effective  at  higher  impact 
energies  when  the  propensity  for  crack  propagation  is  also 
higher.  Comparing  Examples  1 and  4 in  Figure  4.9  one  can  see 
this  phenomenon.  In  Example  4 the  Ingsact  energy  was  very  low 
so  that  stitching  was  not  necessary.  However  in  Example  1 
Che  effectiveness  of  stitch  density  could  be  inferred. 

Prom  Che  above  discussion  we  can  arrive  at  the 
following  conclusions: 

1.  Static  simulations  of  delaminated  stitched  beam 
provide  an  estimate  of  Che  apparent  fracture 
toughness  of  the  stitched  laminates.  The  stitch 
density  significantly  affects  the  increase  in 
apparent  fracture  toughness.  The  percentage 
increase  is  more  for  laminates  with  lower  inherent 
fracture  toughness. 

2.  The  impact  force  at  which  the  delamination  begins 
Co  grow  is  not  dependent  on  stitching.  However 
after  the  delaminacion  growth  is  initiated  stitches 
come  into  play,  and  Che  extent  of  delaminacion 
growth  depends  on  the  stitching  parameters.  In 
general  the  extent  of  crack  propagation  at  the  end 


small  Cor  highsr  stitch 


densities. 

Stitching  is  effective  when  the  impact  energies  are 
higher  and  the  propensity  for  delamtnation  growth 
is  also  higher. 

Stitching  is  more  effective  when  the  delamination 
is  at  the  center. 
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F*impact  force 
q=deflect  at  center  point 
G==energy  release  rate 


Figure  4.6  Flow 
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Center  deflection  q {m] 


4.1  Contact  force  F vs.  center 
(G„5-S30J/!h') 


deflection  q 
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Figure  4.B  Impact 
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□IniClal  crack  length  oC  a apeciaen 

Brinal  crack  length  of  an  unatltched  specimen 

□rinal  crack  length  of  a specimen  with  stitching  (16  ssl) 

C9Flnal  crack  length  of  a specimen  with  stitching  (64  ssl) 


: growth  in  different  examples 
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Crac)c  growth  in  different  examples 
(Gi,c"300J/m*) 


CHAPTER  5 

CONCLUSIONS  AND  FUTURE  WORK 
5.1  SUM.ary 

In  this  dissertation  analytical  and  finite  element 

thickness  stitching  on  the  mechanical  behavior  of 
delaminated  composite  laminates.  In  particular  the  buckling 
and  post-buckling  behavior  and  low-velocity  impact  response 
of  delaminated  composites  are  studied.  These  two  loading 
modes  were  selected  as  they  are  very  detrimental 
to  effective  functioning  of  composites  containing 

The  analytical  models' use  the  Baleigh-Ritz  method  in 
which  stitches  are  modeled  as  smeared  foundation  springs. 
The  finite  element  models  have  the  capability  of  modeling 
individual  stitches.  The  low-velocity  impact  response 


1.  From  the  analytical  models  a single  non-dimensional 
parameter  K'  was  identified: 

. A,E,  NB 
~ Vi  E 

It  may  be  noted  that  K'  is  actually  the  ratio  of  the 
stitch  stiffness  and  Che  equivalent  Young's  modulus  of 
the  beam.  Thus  one  needs  to  use  high  density  stitches 
with  stiffer  yarn  material  if  the  flexural  modulus  of 
the  laminate  is  also  higher.  It  is  found  that  stitching 
increases  the  bucXling  load,  and  reduces  the  energy 
release  cate  at  the  crack-tip. 

2.  The  FE  models  also  support  the  conclusions  of  the 
analytical  models.  The  strains  in  the  stitches  depend  on 
the  axial  stiffness  of  Che  stitches,  and  they  determine 
if  the  stitch  will  break  during  the  post-buckling 
process  and  thus  reduce  the  effectiveness  of  stitching. 

3.  From  the  static  simulations  of  a delaminated  beam  it  is 
found  that  Che  percentage  increase  in  apparent  fracture 
toughness  is  more  for  laminates  with  lower  inherent 
fracture  toughness.  The  impact  force  at  which  the 
delamination  begins  to  grow  is  not  dependent  on 
stitching.  However  Che  extent  of  delamination  growth 
depends  on  the  stitching  parameters.  Stitching  Is 


effective  when  the  inpact  energies  are  higher  and  when 
the  delamination  is  at  the  center. 


5.2  Suggestions  for  Future  Wor)c 

The  present  research  assumed  a single  central 
delanination.  The  models,  especially  the  FB  models,  can  be 
extended  to  beams  with  multiple  delamlnatlons.  Further,  FE 
models  can  handle  laminated  plates  with  stitching  and  other 
forms  of  translaminar  reinforcements.  Since  it  will  not  be 
possible  to  model  individual  stitches  in  larger  plates, 
equivalent  stitch  parameters  such  as  foundation  spring 
constants  may  have  to  be  used.  More  complex  constitutive 
relations  based  on  experimental  observations  can  be  used  to 
model  the  stitch  yarn  or  other  translaminar  reinforcements. 

Since  stitching  results  in  increased  cost  and  weight, 
optimization  studies  can  be  undertaken  to  determine  the 
effective  stitch  materials  and  stitch  pattern  in  the  design 
of  laminated  composite  structures. 


APPENDIX 

DERIVATION  OF  FORCE-DISPLACEMENT  RELATIONS 
IN  THE  BEAM  ELEMENT 


Fxi,  Fzi  and  nomencs  Ci  described  in  Chapter  4 are  given; 


I-. 


Consider  Part  <1>  of  the  beam  shown  in  Figure  A.l.  In 

moment,  po  is  shear  force  provided  by  stitches.  The 
relations  between  external  forces  and  internal  force 
resultants  at  left  side  of  a beam  and  right  side  of  a beam 
are  given  as  follows; 


= -M, 
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= P, 

P\>  * ''i 

C,=M, 


{AD 


Kow  we  consider  Che  internal  force  P in  the  axial  direction 
shown  in  Figure  A. 2. 


The  differencial  equations  of  equilibrium  a 


Here  Af  is  area  of  cross  section  of  part  <1>,  Ei  is  young's 
modulus  of  part  <1>.  Siibstituting  boundary  condition  at 
x«0,  u«ui,  at  x*a.  u»u;  into  Equations  A. 3 and  A.  4 and 
consider  the  equations  in  Equations  A.l  and  A.2r  we  have 
following  expressions  of  f»i  and 


(^5) 


We  consider  the  segment  shown  in  Figure  A. 3 for  deriving 
the  expressions  for  F,,  and  Ci  (1»1,2): 


The  moment 
arbitrary  cross 


p. 


Figure  A. 3 Free  body  diagram 
in  an  arbitrary  section 


and  force  resultants  ff  and  V In  an 
section  can  be  expressed  as: 


- =W  = jW,  +t',x-, 


T:  is  moment  of  ii 


modulus  of 
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(/(9) 


Solving  Equations  A. 7 and  A.S  and  using  the  boundary 
conditions  in  Equations  A. 9 and  A. 10,  we  obtain  relations 
between  V,,  ft,  and  w, , S,.  Then  Pn  and  Ci  can  be  expressed  in 


--V,  -w.  - 


Free  body  diagram  of  portion  <2>  is  shown  in  Figure 


Figure  A. 4 Free  body  diagram  in  part  <2> 

The  procedure  for  deriving  the  expressions  for 
and  Cl  (i=3,4)  is  similar  to  that  one  in  part  <1>,  therefore 
the  derivation  is  omitted  here. 


W.15) 

U.I6) 


_L(_£2*!£l+„  yi° 

d.a  24£,/,  ' 2 
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ASF.J,A,  A 


'2-4.' *'^’'4^,  ^ ^ * 

, a EJz  K P.*i« 


'AA^ 


Ai  is  area  of  cross  section,  Bj  is  Young's  modulus  and  Gi  is 
shear  modulus  in  part  <2>  and  4?'  is  given  by: 


The  free  body  diagram  of  part  <3>  is  shown  in  Figure  A. 5: 


Free  body  diagram  in  part  <3> 


F„,  and  Cl  11-5.6)  . 


A,L  2A,  A,L  2A, 


A,L  2A;  A,L 
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